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ABSTRACT 

Solar active region (AR) 12192 of October 2014 hosts the largest sunspot group in 24 years. It is the most 
prolific fiaring site of Cycle 24 so far, but surprisingly produced no coronal mass ejection (CME) from the 
core region during its disk passage. Here, we study the magnetic conditions that prevented eruption and the 
consequences that ensued. We find AR 12192 to be “big but mild”; its core region exhibits weaker non¬ 
potentiality, stronger overlying field, and smaller fiare-related field changes compared to two other major flare- 
CME-productive ARs (11429 and 11158). These differences are present in the intensive-type indices (e.g., 
means) but generally not the extensive ones (e.g., totals). AR 12192’s large amount of magnetic free energy 
does not translate into CME productivity. The unexpected behavior suggests that AR eruptiveness is limited by 
some relative measure of magnetic non-potentiality over the restriction of background field, and that confined 
fiares may leave weaker photospheric and coronal imprints compared to their eruptive counterparts. 

Subject headings: Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: magnetic topology — Sun: 
photosphere — Sun: surface magnetism 


1. INTRODUCTION 

The great solar active region (AR) 12192 of October 2014 
harbors the largest sunspot group in 24 years. It is so far the 
most intensely flaring region of Cycle 24, producing a total 
of six GOES X-class flares and a multitude of smaller ones. 
Statistically, both flares and coronal mass ejections (CMEs) 
tend to occur in these intense events. A survey of 1996- 
2005 indicates that 75 of 90 X-class flares are associated with 
CMEs; for those above X3 the rate is 23 in 24 (Yashiro et al. 
2006). It is therefore surprising that no CME was detected 
from the core region of AR 12192 during its disk passage. 
Only one on-disk jet-like CME originated from the AR pe¬ 
riphery; two others have been reported to erupt from over 
the east limb (West & Seaton 2015). The unexpected behav¬ 
ior quickly raised interest from the community (e.g., RHESSI 
Science Nugget #239; Thalmann et al. 2015). 

We dub all flares without a CME “confined”: they either 
produce no eruption, or eruptions that fail to escape the Sun 
(e.g., Ji et al. 2003). Confined flares are rare for the more 
energetic events; only a dozen or so confined X-class flares 
have been analyzed (Schmahl et al. 1990; Eeynman & Hund- 
hausen 1994; Gaizauskas et al. 1998; Green et al. 2002; Wang 
& Zhang 2007; Liu 2008; Cheng et al. 2011 ; Chen et al. 2013; 
Liu et al. 2014a). On the other hand, the lack of an associated 
CME is common among weaker ones, accounting for 40% of 
M-class flares (Andrews 2003) and a majority for C-class and 
below. ARs with large eruptive flares often produce smaller, 
confined events too. 

The magnetic cause of confined flares has been studied 
along with the eruption mechanism. Comparative case stud¬ 
ies have largely probed two aspects. One focus is AR non¬ 
potentiality as the source for eruption, e.g., magnetic he- 
licity (Nindos & Andrews 2004) and kink instability (Guo 
et al. 2010). The other focus is the constraining effect of 
the background field, e.g., its decay with height (Liu 2008; 
Guo et al. 2010; Cheng et al. 2011; Nindos et al. 2012) and 


its strength (Wang & Zhang 2007; Liu 2008). Numerical 
experiments seem to suggest that both non-potentiality and 
the background field contribute: by fixing one and adjusting 
the other, confined events can transition to eruptions (Amari 
et al. 2003; Torok & Kliem 2005). Statistical studies have 
shown good correlation between CME occurrence and mag¬ 
netic twist, electric current, and global free energy proxies 
(Ealconer et al. 2002, 2006). Recent work suggests an “upper 
limit” on free energy, where major flares and CMEs preferen¬ 
tially occur (Ealconer et al. 2009; Moore et al. 2012). 

The magnetic consequence of confined flares is little ex¬ 
plored in comparison. A CME bodily removes twisted mag¬ 
netic structure, resulting in decrease of magnetic helicity. The 
flare-related, stepwise change in the photospheric field (e.g.. 
Sun et al. 2012; Wang et al. 2012) have been interpreted to 
be a record of the Lorentz force impulse that drives the ejecta 
(Eisher et al. 2012). Without a CME, do confined flares yield 
smaller magnetic field changes compared to their eruptive 
counterparts? 

Closely watched by multiple observatories, AR 12192 is a 
wonderful test case with all its peculiar behaviors. Here, we 
utilize photospheric vector magnetic data from the Helioseis- 
mic and Magnetic Imager (HMI; Schou et al. 2012) aboard 
the Solar Dynamics Observatory {SDO) to probe the cause 
and consequence of the largest confined flare (X3.1) from AR 
12192. To this end, we select two additional ARs (11429 and 
11158) with major eruptive flares and compare their pre- and 
post-explosion magnetic conditions. Several distinctive fea¬ 
tures of AR 12192 immediately stand out. We discuss the 
implications of our observation. 

2. OBSERVATION AND ANALYSIS 

AR 12192 was well developed when it rotated into view on 
2014 October 17. Subsequent photospheric evolution mainly 
involves fast sunspot separation at two locations (Eigure 1(a)). 
Significant flux emergence commenced at the northern site 
when the AR was near central meridian. The sunspot area 
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Figure 1. Overview of AR 12192. (a) (b) HMI continuum intensity before two X-class flares with scattered light removed. Arrows in (a) denote two locations 
where signiflcant sunspot separation took place. AIA 1600 A ribbons for the X3.1 flare are overplotted in (b). (c) Negative composite AIA image during the X3.1 
flare. Cyan, yellow, and magenta show 131, 171 and 335 A passbands, respectively. Dotted curve outlines the hot 131 A loops connecting the main flaring site to 
the southwest. Boxes in (b) (c) deflne the extent of Figures 2(a). (d) GOES 1-8 A flux (black) and unsigned magnetic flux ^ (green symbols). The AR passed 
central meridian on October 23. CME production of each major flare is marked along the top. Horizontal error bars of indicate 1-day window; vertical error 
bars indicate daily standard deviation. 


(HMI continuum intensity below 0.9 quiet Sun value) reached 
a maximum of ^4300 /iHem. The unsigned flux <l> (sum of all 
pixels where held strength B > 200 G) reached ~2x 10^^ Mx 
on October 27 (Figure 1(d)), an order of magnitude greater 
than the typical total net flux in the polar region during ac¬ 
tivity minimum (Sun et al. 2015). The typical ratio of the 
total signed to unsigned flux is 0.08, suggesting a largely 
closed-fleld environment. All of the conflned flares above 
M3 took place in the core region and showed double chro¬ 
mospheric flare ribbons (Figure 1(b)). The overlying loops to 
the southwest appear to be directly involved in many events 
(Figure 1(c)). 

We aim to compare the magnetic condition of AR 12192 
around SOL2014-10-24T21:41 (X3.1) with several other ARs 
of Cycle 24 that produced major eruptions. We search the 
GOES flare list between 2010 and 2014 using the following 
criteria: their peak intensity is greater than X2; they occur 
between E40 and W40; they produce wide CMEs (halo or 
width greater than 60°); and they show clear, extended dou¬ 
ble ribbons. The last criterion reduces the effect of complex 
magnetic topology; it rules out two candidates with fan-spine 
structure and circular-shaped ribbons (e.g.. Sun et al. 2013). 
Finally, the two regions selected are AR 11429 at SOL2012- 
03-07T00:24 (X5.4) and AR 11158 at SOL2011-02-15T01:56 
(X2.2). The former is the second-most flare-productive AR of 
Cycle 24 (e.g., Wang et al. 2014; Liu et al. 2014b). The latter 
produced the first X-class flare of Cycle 24 and is well stud¬ 
ied (e.g.. Sun et al. 2012). All flares above M3 from these 
two ARs were associated with CMEs. Both ARs also pro¬ 
duced multiple weaker, conflned flares. The opposite-polarity 
sunspots in both ARs underwent strong shearing, in contrast 
to AR 12192. 

HMI generates full disk photospheric vector magnetograms 
with 0.5" plate scale at 12-min cadence. Each magnetogram 
is made from flltergrams taken within a ^22.5 min averag¬ 
ing window, which has a uniform weight for the center and 


less contribution from the edges (Hoeksema et al. 2014). ARs 
are automatically identified and extracted; de-projected maps 
are provided in cylindrical equal area coordinate (Bobra et al. 
2014). Here, we use five maps prior to the onset of each event 
to represent the pre-flare condition, and another five one hour 
after for the post-flare condition. Considering the wide aver¬ 
aging window, we make sure that the last pre-flare maps are 
centered at least 8 min before the flare onset. All contributing 
flltergrams are thus taken well before the hard X-ray peak to 
avoid possible artifacts from the intense flare emission (Qiu 
& Gary 2003). When possible, we consider only strong held 
regions where B is greater than 200 G; results are presented 
as mean ± standard deviation of the five measurements. 

We study the coronal held by extrapolating a nonlinear 
force-free held (NLEEE; Wiegelmann 2004; Wiegelmann 
et al. 2006) and a potential held (PE; e.g., Sakurai 1989) from 
each vector magnetogram (Sun et al. 2012). Magnetic energy 
is calculated as £" = '^B^SV/Stt, where 6V is the grid vol¬ 
ume. The difference of the NLEEE energy En and the PE en¬ 
ergy Ep indicates the free energy Ef that is available to power 
the explosions. We note that should be taken as an order of 
magnitude estimate as the systematic uncertainty can be large 
(Sun et al. 2012). The decay index, n = -dlnBh/dlnz, charac¬ 
terizes the decrease rate of the horizontal held Bh with height 
z (e.g., Kliem & Torok 2006). We also compute the squashing 
factor Q (Demoulin et al. 1996; Pariat & Demoulin 2012) to 
highlight the topological boundaries in the modeled held. 

We have designed a “flaring polarity inversion line” (EPIL) 
mask to demarcate the AR core held, where most free en¬ 
ergy resides (Eigures 2(a)-(c)). We first identify the polarity 
inversion line (PIE) pixels from a smoothed vertical held B^ 
map, and dilate them with a circular kernel (radius r = 3.5 
Mm). Then, we isolate flare ribbons from the 1600 A image 
(above 700 DN s“^) taken near the flare peak by the Atmo¬ 
spheric Imaging Assembly (AIA; Lemen et al. 2012) aboard 
SDO and dilate them with a large kernel. The intersection of 
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Table 1 

Comparison of magnetic characteristics of three major active regions. 





AR 12192 

AR 11429 

AR 11158 

Unit 

Type* 



Flare index 

2335 

1295 

592 





Major flares 

15 

7 

3 



u 


Event 

SOL2014-10-24T21:41 

SOL2012-03-07T00:24 

SOL2011-02-15T01:56 





Location 

S21W21 

N18E31 

S20W10 





GOES class 

X3.1 

X5.4 

X2.2 





Duration 

66 

38 

22 

min 




CME 

No 

Halo 

Halo 





Sunspot area 

4002±11 

1490±2 

861±4 

/iHem 

E 


13 

5r! 


16.12±0.08 

4.88±0.04 

2.73±0.04 

10^2 Mx 

E 


<D 

> 

o 

1 

25.98±0.00 

8.00±0.00 

6.31±0.00 

10'^ A 

E 



logR 

5.30±0.01 

5.32±0.01 

4.89±0.01 


E 

<u 

A 


Mask area 

357±4 

231 ±2 

224±10 

/rHem 

E 

o. 

o 


T> 

0.31±0.03 

0.45±0.03 

0.28±0.04 

1022 Mx 

E 

o 

A 


/ 

l.lOiO.OO 

1.35±0.00 

1.12±0.00 

10‘3 A 

E 



B (rms) 

453±3 

827±4 

678±14 

G 

I 


Ki 

Ph 

J (rms) 

16.1±0.3 

30.7±0.6 

27.7±0.4 

mA 

I 



Shear 

50.2±0.3 

57.9±0.3 

60.1±0.7 

degree 

I 



i«i 

0.05±0.00 

0.19±0.00 

0.25±0.01 

Mm-i 

I 



Current helicity 

0.72±0.04 

6.74±0.06 

5.28±0.20 

G^m-i 

E 


a 

Bum 

220±8 

61±7 

42±0 

G 

I 



Bhm/Bh{2) 

0.35±0.04 

0.06±0.00 

0.05±0.00 


I 

o 

Vh 

o 

o 

> 

O 

Critical height 

77±1 

34±0 

42±1 

Mm 

I 

1 


Ep 

152.8±0.2 

20.9±0.1 

8.8±0.0 

10^2 erg 

E 


<D 

Ef 

4.5±0.0 

10.6±0.0 

2.5±0.0 

10^^ erg 

E 



Ef/Ep 

0.03±0.00 

0.51±0.02 

0.28±0.01 


I 



AEf 

-0.90 

-1.58 

-0.26 

o 

CTQ 

E 



A(J2BhSA) 

+l 

+14 

+8 

10^0 Mx 

E 

tkO 



+11 

+200 

+129 

G 

I 

U 


AFz 

+0.2 

+11.1 

+4.7 

10^^ dyne 

E 



Topology change 

Small 

Large 

Large 


I 


Notes. 

* Indices are classified as extensive (E) or intensive (I); see Section 3.1. Indices in bold are arbitrarily selected as examples for each category. 

^Flare index is defined as +where Mx indicates the GOES magnitude of each X-class flare, etc. Major flares 

include those above M3, between E70 and W70. 

^Sunspot area is computed from HMI intensity, including both umbrae and penumbrae. R measures the total unsigned flux within 15 Mm of 
high-gradient PIL (Schrijver 2007), here with Bz instead of line-of-sight maps. Mean shear is the mean angle between the observed and the 
modeled PE on the photosphere; mean torsional parameter a is calculated as current helicity is approximated by | J2^zJz\ 

(Bobra et al. 2014). 

II The overlying field refers to directly above the FPIL in the PE model. 5/j(42) indicates mean B^ at 42 ± 1 Mm, typical height of eruption 
onset (Liu 2008). BhiAl)/B^il) is the mean ratio of B^ at 42 ± 1 and 2 ± 1 Mm (cf. Wang & Zhang 2007). The critical height is where the B^ 
decay index n reaches 1.5 so the torus instability may set in (Kliem & Torok 2006). 

§The change of the surface integral J^Bh^A and the mean (Bh) consider the FPIL region only, where SA is the pixel area. The change of 
“Lorentz force” Fz refers to the change of ^{B^-B\)5A/{^tt) within FPIL (Fisher et al. 2012). Topological change is assessed qualitatively 
based on Q and coronal field connectivity (Figures 3). 


the dilated PIL and flare ribbons constitutes our FPIL mask. It 
resembles the mask in Schrijver (2007), but includes only the 
part directly involved in a particular flare. Our conclusions 
are not affected if we adjust the mask width (2r) between 5 
and 15 Mm. 

In the largest sunspot umbrae, the HMI inversion module 
sometimes returns unreasonable held values with high formal 
errors. For example, a small patch of abnormally weak 
appeared at the center of the negative sunspot in AR 12192 
(Figure 2(a)). The reason for these “bad pixels” is not fully 
understood; it appears to be the combined effect of low inten¬ 
sity, extremes in the SDO orbital velocity, and limitations of 


the inversion technique. To estimate the adverse effect on our 
analysis, we identify these pixels by setting empirical thresh¬ 
olds on the formal errors and smoothly interpolate over them 
using the data nearby. The difference between the original 
and the interpolated data is 4% for 1% for modeled En and 
Ep, and 8% for (median in time). None affects our conclu¬ 
sions. 

3. RESULTS 

3.1. Pre-Elare Conditions 

We summarize the pre-flare magnetic conditions and the 
flare-related changes for the three ARs in Table 1. Various 
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Figure 2. Comparison of pre-flare magnetic conditions of AR 12192, 11429, and 11158 prior to their respectively largest flare, (a)-(c) maps of the core region. 
They are 164, 120, and 76 Mm respectively in width. The yellow shaded regions denote our FPIL masks, (d)-(f) Maps of the vertically integrated J over the 
lower 11 Mm in the NLFFF model, (g)-(i) Height profile of (black) and decay index n (green) above the FPIL in the PF model. A total of 1-2 x 10^ profiles 
are evaluated for each AR; outliers are removed by using a K-mean algorithm. Lines show the median; shaded bands indicate Icr spread. Horizontal dotted line 
indicates the critical value n = 1.5. In (i), the kink (cf. Nindos et al. 2012) and the larger spread of n are due to the quadrupolar nature of AR 11158 (only the 
central bipole is shown): magnetic connectivity changes rapidly at 10-40 Mm. 


indices shown to be useful indicators for flare and CME ac¬ 
tivity are computed and can be classified as either extensive 
or intensive, following Welsch et al. (2009). Extensive in¬ 
dices generally scale with AR size (e.g., totals) while inten¬ 
sive indices do not (e.g., means). We have the following ob¬ 
servations for AR 12192 regarding the pre-fiare conditions, in 
comparison with ARs 11429 and 11158. 

1. Its global, extensive-type indices are significantly greater. 
These include sunspot area, total magnetic fiux electric 
current /, magnetic energy Ep and and EPIL mask size. 

2. Its extensive-type indices in the core field are comparable 
to the other two ARs. These include the R parameter (a 
free energy proxy, Schrijver 2007), ^ and I within the EPIL 
mask, and free energy Ef. 

3. Its intensive-type indices, particularly those regarding AR 
non-potentiality in the core field, are significantly weaker. 
These include rms field B, rms electric current density /, 
mean shear angle, mean torsional parameter a within the 
EPIL mask, and relative free energy Ef/Ep. The net cur¬ 
rent helicity of the extensive type is small too. This is 
nicely illustrated by the vertically integrated J maps from 
the NLEEE model (Eigures 2(d)-(f)). 

4. Its background field straddling the EPIL is significantly 
stronger. In the PE model, Bk of AR 12192 is stronger in 
low corona (z ^ 42 Mm); the relative strength with respect 
to the near-surface (z ~ 2 Mm) value is higher too. Below 
120 Mm, Bh decreases much slower with height, leading to 
a lower decay index n (Eigures 2(g)-(i)); n does not reach 
1.5 until a large altitude, so the torus instability is less likely 


to set in. We obtain the same conclusions using the NLEEE 
model. 

3.2. Elare-Related Changes 

We evaluate the fiare-related changes using selected indices 
from the last frame before the fiare onset and the frame 1 hr 
after. The actual change near the EPIL is permanent and likely 
occurs on a time scale of minutes (Sudol & Harvey 2005). AR 
12192 is distinctive in the following aspects. 

1. Its photospheric field change is significantly weaker. In the 
two eruptive cases, Bh increased by hundreds of gauss in 
the AR core; such change is not present in AR 12192 (Eig¬ 
ures 3(a)-(c)). To assess the significance, we take the differ¬ 
ence between pairs of Bh maps and compute the rms (a) in 
the EPIL. The pairs are taken both before or after the fiare, 
thus represent a baseline from secular evolution. We find 
that the changes of mean Bh in the EPIL of ARs 11429 and 
11158 reach 2.9cr and 2.2cr respectively (over 5cr in the cen¬ 
tral part), while for AR 12192 the change is 0.2cr. Similarly, 
the change of E^ oc ^(BI-B^), which possibly correlates 
with the ejecta momentum in eruptive fiares (Eisher et al. 
2012), is much smaller in AR 12192. 

2. Its inferred coronal field change is smaller. We compute 
the squashing factor Q from the NLEEE model on a verti¬ 
cal cross section (Eigures 3(d)-(f)). In all three ARs, the 
pre-fiare high-2 patterns indicate the existence of twisted 
magnetic structure (Eigures 3(g)-(i)). In the two eruptive 
cases, such structure largely disappears after eruption, indi¬ 
cating fundamental topological change. The change is less 
significant for AR 12192. 
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Figure 3. Magnetic changes over 1 hr. (a)-(c) Differenced photospheric Bh. Contours show the FPIL masks, (d)-(f) Logarithm squashing factor Q on a vertical 
cut in the NLFFF model, before (left) and after (right) each event. The location of the cut is marked by a short double line in (a)-(c). (g)-(i) Selective field lines 
demonstrating the connectivity changes before (left) and after (right). The horizontal and vertical axes in (a)-(f) have the same scale; height in (g)-(i) is stretched 
for clarity. 


We note that the changes of other indices in Table 1 may 
exhibit a variety of behaviors^. For example, changes very 
little, while current helicity increases significantly in all three 
ARs. We defer investigation of these behaviors to future work. 

4. INTERPRETATION AND DISCUSSION 

AR 12192 exhibits weaker non-potentiality and stronger 
overlying field in the core region, consistent with previously 
studied confined cases. However, the region’s fiare produc¬ 
tivity is extraordinary; most of its extensive properties are 
comparable to or greater than the eruptive ARs 11429 and 
11158. The estimated free energy is enough to power mul¬ 
tiple X-class fiares. We thus argue against global, absolute 
measures as the physical controlling factor of AR eruptive¬ 
ness, and suggest using instead some relative measure that 
quantifies the ratio of magnetic non-potentiality to the restric¬ 
tion of background field. The possibly different controlling 
factors of fiare and CME productivity provide a viable expla¬ 
nation to the strange behavior of AR 12192. 

The exact formula of such a relative measure is unknown 
but may be explored through surveys. The prediction capa¬ 
bility, of course, remains probabilistic. If the measure is rel¬ 
atively low, any fiare from the AR is likely confined. If it is 
high, major fiares are likely eruptive; smaller fiares, however, 
can still be confined. 

HMI observations have confirmed sudden and permanent 
photospheric field change as a common feature in large fiares 
(Sun et al. 2012; Wang et al. 2012). Bh generally strengthens 
along the FPIL, consistent with the magnetic “implosion” sce¬ 
nario where coronal loops contract in response to the reduc¬ 
tion of magnetic pressure from energy release (Hudson 2000). 

Along this line, Fisher et al. (2012) argue that one can use 
the photospheric integral of Maxwell stress tensor to estimate 
the total Lorentz force in a carefully selected volume. The 

^See http://purl.stanford.edu/ss540kq5576 for an extended table. 


temporal integral of its change, which results from the change 
of the photospheric field, can represent the Lorentz force im¬ 
pulse that provides the CME momentum. As the ejecta mo¬ 
mentum is small in the confined X3.1 fiare, the weak change 
of Bh seems to support the argument. 

It is unclear whether all confined fiares exhibit similar small 
changes, although such is consistent with theoretical expec¬ 
tations. In eruptive fiares, the CME bodily removes mag¬ 
netic helicity by ejecting twisted fiux ropes into interplane¬ 
tary space, resulting in less-sheared post-fiare loops (Priest & 
Forbes 2002). In confined events, however, helicity is largely 
conserved. The post-fiare topological complexity is expected 
to differ less, which is what we find. 

The EUV and X-ray observations of AR 12192 show some 
unorthodox features. For example, they have relatively long 
X-ray duration: the X3.1 confined fiare lasted 66 minutes, 
much longer than the two eruptive fiares in ARs 11429 and 
11158 (Table 1); the median duration of all confined fiares 
above M3 is 53 minutes. An estimate of the non-thermal elec¬ 
tron energy for an earlier, confined XI fiare yields ^10^^ erg, 
significantly larger than that of a typical, eruptive XI fiare 
(Thalmann et al. 2015). In addition, the large overlying loops 
that connect to the southwest plage region appear to be con¬ 
tinuously energized (Figure 1(c); see also, Liu et al. 2014a). 
These features suggest a different energy partition in confined 
events. A detailed analysis, similar to those performed on 
eruptive fiares (e.g., Emslie et al. 2004; Eeng et al. 2013), will 
be worthwhile and may shed light on the triggering mecha¬ 
nism. 

We note that our proposed relative measure is different from 
the existing, statistical CME predictors that are extensive in 
nature (e.g., proxies for E’j, Ealconer et al. 2006). There is no 
direct conflict, because AR 12192 is clearly a statistical out¬ 
lier. In fact, any new index must statistically correlate with es¬ 
tablished CME predictors to have any predictive power. Since 
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the Ef proxies also correlate with <l> (Falconer et al. 2002), 
our relative measure should not be scale-free. 

The HMI database has accumulated 4.5 years of vector field 
data, totaling over 4000 strong-field regions and 1.7 million 
records. A comprehensive statistical study is now possible 
and has been performed for fiare likelihood (Bobra & Cou- 
vidat 2015). A similar survey can help clarify whether AR 
12192 is a representative case for confined fiares and can pro¬ 
vide useful insights to CME forecasting. 

We are grateful to Thomas Wiegelmann for the extrapo¬ 
lation code, Seiji Yashiro for the CME-fiare statistics, and 
Anna Malanushenko for the loop-fitting software. This 
work is supported by NASA contract NAS5-02139, awards 
NNX11AJ65G, NNX13AK39G, NSF awards AGS-1321474, 
and AGS-1249150. The SDO data are courtesy of NASA, the 
SDO/HMl and AIA science teams. Magnetic field lines are 
visualized using VAPOR. 
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